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Abstract. This paper presents a counterexample for the approximation 
algorithm proposed by Durocher and Mehrabi [T| for the general problem 
of finding a rectangular partition of a rectilinear polygon with minimum 
stabbing number. 


1 Introduction 

Given a rectilinear polygon P and a rectangular partition R of P, a segment 
is said to be rectilinear relative to P if it is parallel to one of P’s sides. Let 
s be a maximal rectilinear line segment inside P. The stabbing number of s 
relative to R is defined as the number of rectangles of R that s intersects. The 
stabbing number of R is the largest stabbing number of a maximal rectilinear 
line segment inside P. The Minimum Stabbing Rectangular Partition Problem 
(MSRPP) consists in finding a rectangular partition P of P having the smallest 
possible stabbing number. Figure [I] illustrates these definitions. 

Variants of the problem arise from restricting the set of rectangular partitions 
that are considered to be valid. One of these variants is called the conforming 
case , in which every edge in the solution must be maximal, i.e., both of its end¬ 
points must touch the border of the polygon. For this problem, in [Tj , Durocher 
et al. propose an integer programming model for the conforming case where there 
are exactly two edges (that can be in the solution) having each reflex vertex as 
endpoint. Thus, there are also precisely two variables associated to each reflex 
vertex. 

In pQ a 2-approximation algorithm is presented for the conforming case of 
partitions of rectilinear polygons with minimum stabbing number. That approx¬ 
imation algorithm is based in a rounding of the variables. In the Conclusion 
section of the article, it is stated that the algorithm could be extended for the 
general case using a formulation described informally and the same rounding 
rules used in the conforming case. 

In this paper we show that the algorithm as described in T] cannot give a 
2-approximation for the general case of the (MSRPP). This is done by means 
of a counterexample to the referred algorithm. 
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Fig. 1: A rectilinear polygon with a rectangular partition of stabbing number 4. The 
dashed lines represent maximal rectilinear line segments inside the polygon. Segment 
r has stabbing number 4 while segment s has stabbing number 3. 


2 IP Models 

The MSRPP can be modelled via integer programming in a number of different 
ways. In this section we present two such models for the general case of MSRPP 
in an attempt to formalize the description given in [T|. But first, we need some 
definitions. 

Let P be a rectilinear polygon, input of the MSRPP. Define as V p the set 
of reflex vertices of P , i.e., those having internal angles equal to 37t/2. Let 
V p be the set of vertices of P that are not reflex. Denote by grid(P), the set of 
all maximal rectilinear line segments in the interior of P having a vertex in V p 
as one of its endpoints. Let V s p be the set of points in the intersection of two 
segments in grid(P). We refer to these points as Steiner Vertices. The points 
that are not in V p or V p and are in the intersection of a segment in grid(P) 
and the border of P compose the set V p . Denote by V p the set resulting from 
the union of all the point sets defined before, i.e., V p = V p U V p U V p U V p . 

Define E p as the set of line segments in the border of P having only two 
points in V p which are its extremities. Any fragment of a segment in grid(P) 
containing exactly two vertices in V p is called an internal edge. The set of all 
internal edges is E p and the set of all edges in P is E p = E p U E p . A subset 
E' p of E p defines a knee in a vertex u € V p U V p if exactly two edges in E' p 
have u as an endpoint and these edges are orthogonal. A subset E' p of E p is 
said to define an island in a vertex u £ V p if only one edge of E' p have u as an 
endpoint. At last, if ua and ub are two edges in E p having a common endpoint 
it, we denote the angle between ua and ub by 8{ua , ub). 

































Now, we can formalize the model described in [lj as follows: 


( RPST ) 

z = min k 


(1) 

subject to 

%ua + %ub ^ 1) 

Vu£ Vr A ua, ub £ Ef, 

(2) 


%ua + %ub %uc ^ 0, 

V u G V S P , V ua, ub, uc G Ef 




with 8(ua,ub) = tt/2, 

(3) 


y < k — l, 

V set, 

(4) 


uv^E? 
uv p| s^0 




%uv £ ® 

V uv G Ef , 

(5) 


k e z. 


(6) 


In the model above, we have one binary variable x uv for each internal edge uv 
in P which is set to 1 if and only if the corresponding edge is in the rectangular 
partition of P. Constraints (0 ensure that the solution does not contain a knee in 
a reflex vertex. Inequalities m impose that the solution does not form a knee or 
an island in a Steiner vertex. Inequalities © relate the x variables with variable 
k, which represents the stabbing number of the solution. As a consequence, 
the objective function ((1J) is to minimize k. Finally, © and © are integrality 
restrictions for the variables. Figure [5] shows an instance of the MSRPP (called 
random-20-17) with 62 internal edges and their corresponding variables. 
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Fig. 2: Instance random-20-17 with 62 internal edges and the corresponding variables. 




















































































As stated before, the ( RPST ) model above is not the only model for the 
problem and next we show another way of modelling it. However, to guarantee 
the correctness of the model we must first prove a property of optimal solutions 
for the MSRPP. The following proposition is a generalization of Observation 1 
in ffl. 

Proposition 1. Any rectilinear polygon P has an optimal rectangular partition 
R in which every maximal segment of R has at least one reflex vertex of P as 
an endpoint. 

Proof. Let R be a rectangular partition of a rectilinear polygon P. Let e be a 
maximal segment in R having a and b as its endpoints. Suppose neither a nor 
b are reflex vertices. Since e is maximal and R is a rectangular partition, both 
endpoints of e must lie in segments perpendicular to e. 

Now, since R is a rectangular partition, e define two minimal rectangles (each 
one possibly containing other rectangles) having e as one of its sides, let us denote 
them by r\ and r 2 - There are three cases to consider. 

The first case consists of r± and r 2 been empty rectangles, i.e., neither rq nor 
r- 2 . contain other rectangles. Therefore, the removal of e unite these rectangles, 
composing a single rectangle. Therefore, R\e is still a rectangular partition. It 
is clear that removing a segment cannot increase the stabbing number of the 
solution. Thus, if R is an optimal solution, so is R \ e. 

The second case to consider is when only one of rq or r 2 contains other rect¬ 
angles. Suppose without loss of generality that rq is the one containing other 
rectangles. Now, we can drag e towards rq, shrinking any segment with an end¬ 
point in e, until e meets a reflex vertex or the border of P. In the latter case, e 
is merged to the border of P. It is easy to see that the result of this dragging 
operation is also a rectangular partition besides, the only stabbing segments af¬ 
fected by this operation are the ones parallel to e and their stabbing number 
cannot increase. Therefore, as R is optimal, so must be the new solution. 

At last, we must consider the case where both rq and r 2 contain other rect¬ 
angles. Suppose without loss of generality that the number of segments in rq 
having an endpoint in e (thus, perpendicular to it) is greater or equal than the 
number of segments with these characteristics in r 2 ■ Then, again, we can drag e 
towards rq, shrinking any segment with an endpoint in e, until e meets either a 
segment parallel to e or a reflex vertex or the border of P. If a parallel segment 
is met, e is merged to it and the process is repeated until a reflex vertex or the 
border of P is met. In case the border of P is met, e ceases to exist together with 
the segments in the space between e and the border. Once again, the dragging 
operation results in a rectangular partition of P and the only stabbing segments 
affected by this operation are parallel to e. But, as the number of segments in 
r 1 is greater or equal than the number of segments in r 2 , one can see that the 
stabbing number of the new rectangular partition cannot be greater than that 
of R. 

Ergo, there is always an optimal rectangular partition where every maximal 
segment has at least one reflex vertex of P as an endpoint. □ 


In the next model, given the same definitions as before, we consider the set 
E p of rectilinear segments uv where u G V p and v G V p . Notice that a segment 
of E p can be comprised of several consecutive segments of E p . Hence, we call 
E p the extended edge set. In the formulation below, we have a variable x uv 
for each edge in E p and from Proposition Q] it is easy to notice that this set of 
variables is sufficient to provide optimal rectangular partitions. 


(. RPST2 ) z 
subject to 


min k 


(7) 

^ ^ %ua ^ I? 
uaeE^ 

v«e v p 

(8) 

Xab + X U v ^ 1) 

V ab, uv : ab (~l uv =£■ 0 A 



A ah n uv ^ a, b, u or v 

(9) 

^ ^ Xuv Xab ^ 0, 

0(uv,ab)=TT/2 A 

A b£uv A b^u A b^v 

if a £ V P , b G Vf 

(10) 

^ ^ Xuv k 1, 

uv^EP :uv p| s^0 

V s G L 

(11) 

Xuv £ B 

V uv G E P . 

(12) 

k GZ 


(13) 


In this model, inequalities (|5]) guarantee that the solution does not contain 
a knee in a reflex vertex. Constraints m enforce planarity (two segments of 
the partition can only intersect at their extremes). Constraints (11011 prevent the 
existence of knees and islands in a Steiner vertex. Finally, CCD are the stabbing 
constraints and da and m are integrality constraints. Figure OH shows instance 
random-20-17 with 42 internal edges and the corresponding variables. 

3 The Counterexample 

Before discussing the counterexample, we first present the rounding scheme pro¬ 
posed in [I] for the conforming case. Once the optimum of the linear relaxation 
is computed, the rules for rounding variables in the conforming case are really 
simple: a variable corresponding to a horizontal segment is rounded down to zero 
if its value is smaller than or equal to 0.5 and is rounded up to one if its value 
is greater than 0.5. A variable corresponding to a vertical segment is rounded 
down to zero if its value is smaller than 0.5 and is rounded up to one if its value 
is greater than or equal to 0.5. 

In the Conclusion section of [I], a model for the general (non-conforming) 
case is described informally. From the discussion, apparently such model is equiv¬ 
alent to the ( RPST ) formulation given in Section [2] According to the authors, 
the same rounding rules used in the conforming case provide a 2-approximation 
for the general case. 
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Fig. 3: Instance random-20-17 with its extended edges and corresponding variables. 



The rounding rules do not mention what should be done for Steiner vertices, 
and no guarantee is given that applying them directly in these situations will 
avoid the formation of a knee or an island. In fact, the instance displayed in 
Figure [I] shows that this cannot always be done without sacrificing feasibility. In 
this figure, the optimal values of the variables corresponding to edges incident 
to Steiner vertex v39 (see Figure [3]) after solving the linear relaxation associated 
to instance random-20-17 are given. As only the variable corresponding to one 
vertical edge incident to that vertex has value equal to 0.5 and the other three 
are smaller than 0.5, rounding according to that rule would result in an island 
at ?;37. Therefore, the set of edges obtaining after rounding does not form a 
rectangular partition. 
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Fig. 4: Values of variables corresponding to edges incident to a Steiner vertex after 
solving linear relaxation. 


It is however possible that we misinterpreted the model the authors were 
thinking of (although there is evidence in contrary) and the idea is actually to 
define variables corresponding to all edges having a reflex vertex as one of its 
endpoints. If so, the formulation would look like (RPST 2) model in the previous 
section. In this alternative formulation, rounding the variables using that rule 
does not cause the same problem as before since every variable correspond to an 
edge having a reflex vertex as endpoint. 

Contrary to what happens in the conforming case, however, the reflex vertices 
here have more than two incident edges. Therefore, it is possible that the solution 
of the linear relaxation result in values smaller than 0.5 for all the variables 
corresponding to the edges incident to a certain reflex vertex. Thus, the rounding 
of such solution would result in a partition having a knee in a reflex vertex. 






















The situation described above occurs in practice with instance random-20-17, 
as shown in Figure [5] Consider the edges incident to vertex v5. All the associated 
variables incident to this vertex have value smaller than 0.5. As consequence, 
they will be rounded to zero, resulting in the formation of a knee at v5 and, 
therefore, in an infeasible solution. 
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Fig. 5: Values of variables corresponding to edges incident to a reflex vertex after solving 
linear relaxation. 


4 Conclusion 

From the counterexample presented in Section [31 we conclude that it remains 
open whether a 2-approximation for the MSRPP in the general case exists. It 
is, however, noteworthy that many other contributions are presented in [T] and 
none of them are diminished by this counterexample. 
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Appendix 


File name: random-20-17.rect 
Model: RPST 
Vertex number: 59 
Edge number: 62 

Reading Problem stab 
Problem Statistics 

231 ( 0 spare) rows 

63 ( 0 spare) structural columns 

752 ( 0 spare) non-zero elements 

Global Statistics 

63 entities 0 sets 0 set members 

Minimizing MILP stab 
Original problem has: 

231 rows 63 cols 752 elements 63 globals 

Its Obj Value S Ninf Nneg Sum Inf Time 

0 .000000 D10 24.000000 0 

74 2.416667 D 0 0 .000000 0 

Optimal solution found 
*** Search unfinished *** Time: 0 
Number of integer feasible solutions found is 0 
Best bound is 2.416667 

Solution: 


XI 

= 

0.583333 

X2 

= 

0.416667 

X3 

= 

0.095238 

X4 

= 

0.916667 

X5 

= 

0.666667 

X6 

= 

0.333333 

X7 

= 

0.428571 

X8 

= 

0.571429 

X9 

= 

0.904762 

X10 

= 

0.845238 

Xll 

= 

1.000000 

X12 

= 

0.571429 

X13 

= 

0.000000 

X14 

= 

1.000000 

X15 

= 

0.000000 

X16 

= 

1.000000 

X17 

= 

-0.000000 

X18 

= 

0.500000 

X19 

= 

0.238095 

X20 

= 

-0.000000 

X21 

= 

-0.000000 

X22 

= 

0.333333 

X23 

= 

0.428571 

X24 

= 

1.000000 

X25 

= 

0.238095 

X26 

= 

-0.000000 

X27 

= 

0.285714 

X28 

= 

0.130952 

X29 

= 

0.500000 

X30 

= 

0.333333 

X31 

= 

0.333333 

X32 

= 

0.333333 

X33 

= 

0.333333 

X34 

= 

0.333333 

X35 

= 

0.095238 

X36 

= 

0.583333 

X37 

= 

0.333333 

X38 

= 

0.500000 

X39 

= 

0.500000 

X40 

= 

0.297619 

X41 

= 

0.285714 

X42 

= 

0.583333 

X43 

= 

0.285714 

X44 

= 

-0.000000 

X45 

= 

0.333333 

X46 

= 

0.333333 

X47 

= 

-0.000000 

X48 

= 

0.285714 

X49 

= 

0.285714 

X50 

= 

0.285714 

X51 

= 

0.285714 

X52 

= 

0.571429 

X53 

= 

0.285714 

X54 

= 

0.285714 

X55 

= 

0.571429 

X56 

= 

0.000000 

X57 

= 

0.714286 

X58 

= 

1.000000 

X59 

= 

-0.000000 

X60 

= 

0.500000 

X61 

= 

0.285714 

X62 

= 

1.000000 

X63 

= 

3.000000 








File name: random-20-17.rect 
Model: RPST2 



Vertex number: 59 
Edge number: 62 
Reading Problem stab 
Problem Statistics 

336 ( 0 spare) rows 

63 ( 0 spare) structural columns 

996 ( 0 spare) non-zero elements 

Global Statistics 

63 entities 0 sets 0 set members 

Minimizing MILP stab 
Original problem has: 

336 rows 63 cols 996 elements 63 globals 

Crash basis containing 13 structural columns created 

Its Obj Value S Ninf Nneg Sum Inf Time 

0 .000000 D10 24.000000 0 

66 2.413793 D 0 0 .000000 0 

Optimal solution found 
*** Search unfinished *** Time: 0 
Number of integer feasible solutions found is 0 
Best bound is 2.413793 

Solution: 


XI 

= 

0.190476 

X2 

= 

0.380952 

X3 

= 

0.142857 

X4 

= 

-0.000000 

X5 

= 

-0.000000 

X6 

= 

-0.000000 

X7 

= 

-0.000000 

X8 

= 

-0.000000 

X9 

= 

0.285714 

X10 

= 

-0.000000 

Xll 

= 

0.142857 

X12 

= 

0.142857 

X13 

= 

-0.000000 

X14 

= 

0.000000 

X15 

= 

0.523810 

X16 

= 

-0.000000 

X17 

= 

-0.000000 

X18 

= 

0.190476 

X19 

= 

0.238095 

X20 

= 

0.095238 

X21 

= 

0.000000 

X22 

= 

-0.000000 

X23 

= 

0.476190 

X24 

= 

-0.000000 

X25 

= 

0.190476 

X26 

= 

0.380952 

X27 

= 

0.190476 

X28 

= 

0.000000 

X29 

= 

0.238095 

X30 

= 

-0.000000 

X31 

= 

-0.000000 

X32 

= 

0.190476 

X33 

= 

0.523810 

X34 

= 

0.809524 

X35 

= 

-0.000000 

X36 

= 

0.285714 

X37 

= 

0.380952 

X38 

= 

0.000000 

X39 

= 

0.000000 

X40 

= 

0.190476 

X41 

= 

0.619048 

X42 

= 

-0.000000 

X43 

= 

-0.000000 

X44 

= 

0.142857 

X45 

= 

-0.000000 

X46 

= 

-0.000000 

X47 

= 

-0.000000 

X48 

= 

0.380952 

X49 

= 

0.000000 

X50 

= 

0.380952 

X51 

= 

-0.000000 

X52 

= 

0.619048 

X53 

= 

-0.000000 

X54 

= 

-0.000000 

X55 

= 

0.095238 

X56 

= 

0.142857 

X57 

= 

-0.000000 

X58 

= 

-0.000000 

X59 

= 

0.380952 

X60 

= 

-0.000000 

X61 

= 

0.000000 

X62 

= 

0.380952 

X63 

= 

3.000000 









